QQ-onia package: a numerical solution to the Schrödinger radial equation for heavy quarkonium

This paper presents the basics of the QQ-onia package, a software based upon the Numerov 0(h) method which can be used to solve the Schrödinger radial equation using a suitable potential V (r) for the heavy quarkonium system. This package also allows the analysis of relevant properties of those resonances such as the square of the wave functions at the origin, their corresponding derivatives for l 6= 0 states, typical heavy-quark velocities, and mean square radii. Besides, it includes a tool to analyze the spin dependent contributions to the heavy quarkonia spectrum, providing the splitting of nS1−nS0, as well as the nPJ −nP1 energy levels. Finally a simple software implemented in QQ-onia to compute E1 transition rates is presented.